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Abstract 

We perform a fit to the recent XenonIOO data and study its implications for Dark 
Matter scenarios. We find that Inelastic Dark Matter is disfavoured as an explana- 
tion to the DAMA/LIBRA annual modulation signal. Concerning the scalar singlet 
DM model, we find that the XenonIOO data disfavors its constrained limit. We 
study the CMSSM as well as the low scale phenomenological MSSM taking into 
account latest Tevatron and LHC data (1.1 /fb) about sparticles and B s — > /i/i. 
After the EPS 2011 conference, LHC excludes the "Higgs-resonance" region of DM 
freeze-out and XenonIOO disfavors the "well-tempered" bino/higgsino, realized in 
the "focus-point" region of the CMSSM parameter space. The preferred region 
shifts to heavier sparticles, higher fine-tuning, higher tan j3 and the quality of the 
fit deteriorates. 

1 Introduction 

The best motivated candidates for the cold Dark Matter (DM) of the Universe [1] are Weakly 
Interacting Massive Particles (WIMPs) [2]. Indeed, a stable particle with a typical weak inter- 
action cross section of 1 pb and a mass of order 100 GeV naturally produces the observed DM 
thermal relic abundance. Because the evidence for DM is purely gravitational, there are many 
particle physics models for WIMPs. Theoretically most appealing ones among those are the 
DM models based on supersymmetry (SUSY) [3] in which the DM is stable due to the existence 
of discrete Z2 symmetry - the R-parity [4] or equivalently the matter parity [5]. Because the 
matter parity may occur also in the scalar extensions of the standard model (SM) [6], the scalar 
DM models are also well motivated from the underlying unified physics point of view. 



1 



The WIMPs have been searched for in the LEP, Tevatron and LHC experiments as well as 
in the underground DM direct detection experiments that detect DM recoils on nuclei. The 
history of DM direct detection experiments can be characterized with their sensitivity to the 
spin-independent DM-nucleon interaction cross section a^i as follows. 

1. The non-relativistic DM annihilation cross section suggested by the observed cosmo logical 
DM abundance is o ~ l/(T Mpi) w 3 10 _26 cm 3 /s = 10 _36 cm 2 . Such a large value of cxsi 
is clearly excluded unless the DM particles have peculiar properties or their interactions 
are restricted to untested sectors. 

2. DM that interacts with the Z boson: crgj ~ a 2 m 2 N /M^ w 10~ 38 cm 2 . Again, generic 
WIMPs of this type have been excluded. 

3. DM that interacts with the Higgs boson: crgi ~ a 2 m^/Mp M M| 10 _43 cm 2 , having 
assumed that the DM-DM-Higgs coupling is comparable with the weak coupling a. 

4. DM that interacts with the W ± boson and consequently at loop level with nucleons [7], 
cr SI f» a 4 m% / (4n) 2 Mv M Myy 10~ 46 cm 2 , as predicted e.g. by Minimal Dark Matter [8]. 

Therefore, in the absence of new gauge bosons, the present most sensitive DM direct search 
experiments like CDMS [9], EDELWEISS [10] and XenonIOO [11] are testing particle physics 
models in which DM is coupled to the Higgs boson or in which the DM interactions with matter 
occur only at loop level. 

Although the CDMS, EDELWEISS and XenonIOO searches for the DM have not given a 
discovery so far, there is a long-standing claim by DAMA/LIBRA [12] that sees an 8a evidence 
for an annual modulation in the nuclear recoil signal of low-mass DM. Recently this claim 
received some support from the results of CoGeNT experiment [13]. In order to reconcile this 
result with the negative searches from other experiments, un-orthodox DM scenarios such as 
Inelastic Dark Matter (iDM) [14, 15] have been put forward. Testing all those scenarios is of 
utmost importance for understanding the nature of DM. 

Recently the XenonIOO experiment published new data after 101 days of data taking 
that probes values of the spin- independent DM/nucleus cross section down to crgi ~ 7.0 x 
10 _45 cm 2 [16]. They observed 3 signal candidate events with the expected background of 
(1.8 ± 0.6) events. This result leads to the most stringent limit on DM interactions today, and 
further constrains the best motivated DM models. 

In this paper we perform a fit to the new XenonIOO data and apply the results to constrain 
several model dependent as well as model independent scenarios of DM. First, motivated by 
the DAMA/LIBRA anomaly we study whether the iDM, as an explanation to the observed 
annual modulation, can survive the new XenonIOO results. We find that this is not the case 
and the DAMA/LIBRA anomaly needs another explanation. After that we study the simplest 
possible DM model, the scalar singlet model. 

However, most of our effort goes to the studies of SUSY models. We study the Constrained 
Minimal Supersymmetric Standard Model (CMSSM) as well as the low energy phenomenologi- 
cal MSSM (pMSSM). In doing that we apply also constraints on SUSY parameter space coming 
from the LHC at = 7TeV with 1.1/fb of integrated luminosity, new data from DO, CMS 
and LHCb on B s — > and new data from Tevatron on the top quark mass [17, 18, 19]. 

We find that the XenonIOO data stringently constrains the "well-tempered" neutralino sce- 
nario [20], and excludes the corresponding "focus-point" parameter region of the CMSSM [21]. 
Therefore we perform a model independent analysis of that scenario. We derive constraints of 
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parameter 


range 


distribution 


Vq kins" 1 


200 - 240 


Gaussian: 220 ± 10 


v esc kms -1 


498 - 618 


see [22] 




±2a 


see [16] 



Table 1: The parameters ranges for the XENONiOO fits. 

new XenonIOO data on the CMSSM and pMSSM parameter spaces and study the implica- 
tions of those to the expected results in the LHC experiments. We find that, after EPS 2011 
conference, the "Higgs-resonance" mechanism of DM freeze-out is completely excluded. The 
allowed CMSSM parameter space is pushed towards higher sparticle masses, the allowed range 
of tan /3 is constrained from below and from above, and the fine tuning of the parameters has 
become more severe. We show that the DM direct detection experiments and the LHC are 
complementary in constraining the supersymmetric parameter space. 

The paper is organized as follows. In section 2 we perform the fit of new XenonIOO data 
and apply this to iDM in section 3. In section 4 we study the scalar singlet model. In section 5 
we study supersymmetric dark matter. We conclude in section 6. 



2 Fit of XenonIOO data 

We perform a fit of the new data relased by the XenonIOO collaboration [16], taken between 
January and June 2010, which corresponds (before any cut) to an exposure of 101 days x 48 kg. 
XenonIOO reports 3 events in the signal region which are consistent with an estimated back- 
ground of 1.8 ± 0.6 events. 

For fixed values of M DM and a^i, the number of signal events expected in a given direct- 
detection experiment depends on astrophysical, nuclear and experimental parameters which are 
more or less well-measured. We build a global x 2 including these extra parameters as nuisances 
and marginalize over them at the end. 

On the astrophysical side one has to know both the local DM density, p & , and the DM veloc- 
ity distribution. In the standard halo model the latter is fixed to be a Maxwellian distribution 
with a given r.m.s. velocity vo, truncated at certain escape velocity ve- 

dN 9 ; 9 

— = f(v)oce- v2 / v o9(v-v E ). (1) 

The preferred values of the two velocities are vq = 220 kms" 1 and ve = 544 kms" 1 [22]. We 
assume that Vq has a Gaussian uncertainty of ±10 kms -1 , and we allow for a 2er variation of it 
in the fit. We take from [22] the probability distribution of the ve parameter. The local DM 
density p is kept fixed to 0.3 GeV /cm 3 since the effect of its variation is just equivalent to 
an overall rescaling of o~si (the quantity really probed by direct detection experiments is their 
product osiPo). 

In terms of this set of parameters the differential scattering rate per unit detector mass is 

N TT7 — / d A vvf[v)—— , (2) 



dE nT M DM J\$\ >Vn . dE„ 



where is the number of nuclei in the target per unit detector mass, Mdm is the DM mass 
and f m ; n is the minimal velocity which allows a nuclear recoil with energy E nT . ctdm n is the 
nucleus-DM cross section. 
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We will be interested in spin-independent processes. Nuclear physics enters into the de- 
termination of form factors which modulates the cross section according to the momentum 
transfer: 

^dm n _ m N a sl (f p Z + f n (A - Z)f J?2f 

Here m N is the nucleus mass, \i n is the nucleon-DM reduced mass, and ogi is the spin- 
independent nucleon-DM cross section. We fix our choice of F to the standard Helm form 
factor [23]. We furthermore assume equal coupling strength to both protons and neutrons: 

fp^fn = 1- 

The last important source of uncertainty comes from the experiment itself. In XenonIOO 
the energy of a recoiling nucleus is inferred from the number of recorded photo-electrons (called 
Si signal) emitted after the prompt relaxation of the excited nucleus. The relation between SI 
and the recoil energy is given by 

Sl(E m ) = 3.6 PE x E nr x C cS (E nT ). (4) 

The value of the function C e g at recoil energies (nuclear) smaller than 5keV is not measured 
and must be extrapolated. Such extrapolation is crucial to put limits on light dark matter: the 
larger C e s is, the stronger the limits are. We use the recent C e s measurement from [16], based 
on [24], and we allow in our fit 2a excursions from its central value. 

In view of the small number of observed events we perform an event-by-event fit, such that 
we construct the x-square [25] 

XXenonl00( M DM, 0"SI, {Pi}) = X{ Pl }(M DM , <T SI ) + ^ xliPi), (5) 

i 

where x\ is the one associated to the parameters themselves and x\ Pi ] is the %-square for fixed 
values of the parameters Pi (v , v esc and £ e ff)> evaluated in terms of the SI signal along the lines 
of [11]. We thus include both a Poissonian smearing of the signal and a Gaussian smearing 
of the energy measurement. The events are collected in the window 4PE < SI < 30 PE. We 
use the acceptance from [16] and we model the background as a flat distribution in the signal 
region normalized to reproduce the expected number of background events. 

We then marginalize Xxenonioo "with respect to the nuisance parameters Pi defining 

XXenonlO()( M DM, 0Sl) = Him xLnonlOO ( M DM, 0"SI, {Pi}), (6) 

{Pi} 

where the minimum is taken in the range summarized in Table 1. 



3 Inelastic Dark Matter 

The DAMA collaboration reported an 8a evidence for an annual modulation of their nuclear 
recoil signal [12] which is compatible with light dark matter scattering on nuclei. Such signal 
conflicts with other direct detection experiments and (as already remarked by the XenonIOO 
Collaboration [16]) is strongly disfavored by XenonIOO under the assumption of an elastic 
interaction of the dark matter with the Na nuclei of the DAMA detector. We show this in 
Fig. 1 where we plot the region favored by elastic scattering on sodium at DAMA together with 
the XenonIOO exclusions. The DAMA confidence levels are obtained by a simple x 2 fit of the 
first 12 bin of fig. 9 in [12], showing the energy dependence of the seasonal modulation. 
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Figure 1: The 95, 99.7% confidence level contours for 2 d.o.f. for the DAMA modulated signal 
under the assumption of elastic scattering on sodium atoms. We assume different values of 
the sodium quenching factor, 0.2, 0.3, 0.4 and 0.5, from right to left. The blue lines are the 
Xenon 100 exclusion curves at 95% (continuous curve), 99.7% (dashed) confidence level. We 
assume vq = 220 km s -1 and v esc = 544 km s -1 for the DAMA fit. We neglect channeling. 

The strong bounds already from previous experiments, prompted theorists to consider In- 
elastic Dark Matter [14] (iDM), namely the idea that the DM detection process could be 
DM iV — > DM' N where the new state DM' is heavier than the DM state by an amount 
5 = M DM / — Mdm- It was pointed out that the modified kinematics could make the DAMA 
anomaly compatible with experiments, such as CDMS, that use lighter nuclei. Indeed the min- 
imum DM velocity needed to scatter on a nucleus with mass mjv giving recoil energy E nr is [14] 



where fi^ is the reduced mass of the DM/nucleus. For the mechanism to be effective the 
inelasticity mass splitting is assumed to be around 100 keV, which is the order of the typical 
nuclear recoil energy. By assuming progressively smaller tails of the uncertain DM velocity it is 
possible to avoid progressively stronger constraints from experiments with light nuclei. For the 
same reason it is possible to enhance the relative modulation in the signal claimed by DAMA. 

Xenon nuclei have a mass very close to iodine which is, at DAMA, the dominant source 
of recoils for DM particles with mass above 0(10 GeV). This fact, together with the high 
exposure of XenonIOO allows to put strong constraints on the iDM hypothesis. This is shown 
in Fig. 2 where we plot 95% and 99.7% confidence level contours for the parameters favored by 
DAMA in the (5, cxsi) plane, together with the XenonIOO exclusion curves. The Vq, v esc and 
£ e ff parameters have been here fixed to their median values. We find that the iDM hypothesis 
to explain DAMA is disfavoured throughout its parameter space. This conclusion is sound, 
independent of the various uncertainties thanks to the similarity between Xe and I nuclei. This 
can be understood in a simple way [15]. The rate observed by DAMA can be written as 



where B is a time-independent unknown background, S is the mean DM signal and a describe 
the size of the modulation, to =2 June, is where the peak of the modulation observed by DAMA 




(7) 



B + S(l + acosu(t — t )), 
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Figure 2: The 95, 99.7% confidence level contours for 2 d.o.f. for iDM fit to DAM A, together 
with the 95, 99.7% exclusion curves from XENONiOO data (full and dashed respectively). We 
fix the iodine quenching factor to 0.085. Sodium quenching and channeling are irrelevant for 
these values of the mass. 

is located. The size aS of the measured modulation is ~ 0.02 counts kg -1 day -1 in the range 
2keV < E nT < 6keV. Due to the iodine quenching factor, qj f» 0.085, this region translates to 
23keV < E nr < 70keV at Xenon. In the region of overlap between the former range and the 
acceptance window 8.4 — 44.6 keV, the DM scattering rate at Xenon is given by 

-RxenonlOO ~ A X € X S (l/fl + COS u(t ~ t )) ■ (9) 

e and A are respectively the efficiency of the XenonIOO cuts and its acceptance. Together they 
amounts to roughly 0.2. Integrating eq. (9) over the energy range and over the time of Xenon 
data-taking one finds 

iVxcnonioo > 0(50) events, (10) 

where we folded in the exposure and we used a < 1. It is thus clear where the exclusion curves 
shown in Fig. 2 comes from. 



4 Scalar singlet Dark Matter model 

We consider a DM model obtained adding to the Standard Model a Dark Matter real singlet 
scalar field S coupled to the Higgs doublet H as described by the following Lagrangian [26, 27, 
28, 29]: 

2 = ^sm + ^£ - ^S 2 - XS 2 \H\ 2 , (11) 

invariant under the matter parity S — > — S, the discrete gauge symmetry of scalars carrying 
B — L [6] . S is the DM field and its mass is given by M£ M = m 2 + XV 2 with V = 246 GeV. So 
the model has 2 free parameters, Mom and A. Assuming that the relic DM abundance equals 
its cosmologically measured value the model is able of predicting the relation in the plane 
(Mdmj °"si) plotted in Fig. 3. Given that the Higgs boson mass is also presently unknown, we 
plotted such relation for a few values of the Higgs mass: 115 GeV (green, favored by precision 
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Figure 3: Predictions of the scalar singlet model for a few values of the Higgs boson mass: 
115 GeV (green), 140 GeV (yellow), 200 GeV (red), 300 GeV (magenta). The dots are the 
predictions of the constrained model of [6, 30}. 

data), 140 GeV (yellow, compatible with precision data), 200 GeV (red, disfavored by precision 
data); 300 GeV (magenta, strongly disfavored by precision data). 

The very small crsi predicted around Mdm = wifc/2 is due to the Higgs resonance enhance- 
ment of the cosmological DM annihilation rate. This quantity is computed as in [30] , already 
including the 3 body final states whose relevance was emphasized in [31]. 

In absence of a theoretical motivation for having m comparable to the Higgs mass, [6, 30] 
considered the case m — 0, such that the model has one parameter less and is able of predicting 
a point in the plane (Mdm, 0si). Such prediction is also shown in Fig. 3, for the same values of 
the Higgs boson mass. 

We remark two uncertainties not explicit from the plot. First, the XenonIOO exclusion 
bound is plotted assuming p Q = 0.3GeV/cm 3 for the local DM density. This is the canonical 
value routinely adopted in the literature, with a typical associated error bar of ±0.1 GeV/cm 3 . 
Recent computations found a higher central value closer to 0.4GeV/cm 3 [32] that would imply 
stronger bounds on the cross section cr SI . 

Second, the prediction for the conventional spin-independent DM/nucleon cross section is: 

asi = ttM 2 m 4 ' ( } 

where / parameterizes the nucleon matrix element: 

(N\m q qq\N)=f q m N [NN], f= £ f q = I + 1 £ f q . (13) 

q={u,d,s,c,b,t} q={u,d,s} 

The main uncertainty comes from f s . The recent analyses use / = 0.56±0.11 [33], or / = 0.30± 
0.015 [34], in agreement with the lattice results [35] and phenomenological determination [36]. 
Here and in the following we assume the default value in the Micromegas code: / = 0.467 [37]. 
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quantity 


experiment 


Standard Model 


a^lMy) [461 


1184 ± 0007 

V J • _L _1_ v7 -L — 1 — V 7 • V / V J V J 1 


parameter 


m+ [471 


173 1 ± 9 

_I_ 1 • _L — 1 — V J • '.7 


parameter 


m b [48] 


4.19 ±0.12 


parameter 


Odm/j 2 [49] 


0.112 ±0.0056 





5^ [50] 


(2.8±0.8)10- 9 





BR(E d -> X s7 ) [51] 


(3.50 ± 0.17) 10- 4 


(3.15 ±0.23) 10~ 4 


BR(B S -»■ [19] 


(0.9 ±0.6) 10~ 8 


(0.33 ±0.03) 10~ 8 


BR(£ U -> rz/)/SM [52] 


1.25 ±0.40 


1 



Table 2: The data we fit, together with LHC and Xenon 1 00 bounds. 

5 Supersymmetry 

In this section we study the impact of new XenonIOO data on constraining SUSY models. We 
first consider the CMSSM, the most popular SUSY model with an unified scalar mass m , an 
unified gaugino mass Mi/ 2 and an unified trilinear scalar /4-term at the GUT scale. Given that 
XenonIOO adds to many other experimental constraints, we perform a global fit to all relevant 
data as described in the next subsection. Most importantly, we include the recent CMS and 
ATLAS constraints on the CMSSM parameters space that are based on the LHC data with 
1.1/fb presented in july 2011. Our fits extend the previous ones obtained in the similar studies 
in [38, 39, 40, 41, 42, 44, 45] without EPS 2011 data. 

One of the results of our research is that the "well-tempered" neutralino scenario is strin- 
gently constrained by the XenonIOO new result. Motivated by that we perform a generic 
analyses of the "well-tempered" neutralino and show that, quite model independently, this 
scenario is now stringently constrained. 

Finally we relax the CMSSM constraints on the particle spectrum coming from the uni- 
fication relations and study a generic low energy phenomenological MSSM, the pMSSM. We 
identify the generic XenonIOO constraints on pMSSM. 

5.1 Global fit of supersymmetric models 

Global fits of the SUSY models have recently been performed by several groups in the context 
of LHC studies, and our results agree with them. Therefore we focus on the impact of the new 
XenonIOO data in constraining SUSY models. Here, we briefly describe our procedure, and 
the differences with respect to previous approaches. 

We perform a random scan of the parameter space and calculate the sparticle spectrum as 
well as DM relic abundance using the Micromegas public code [37]. We vary the relevant input 
SM parameters (m t , a^, m^) within 5a experimental errors assuming a Gaussian distribution. 
The CMSSM parameters (m , Mu2, A , tan/3, sign(^)) are generated randomly in the ranges 
(m ,M 1/2 ) ~ (0,4000) GeV, A ~ (-3m ,3m ), tan/3 ~ (1,60) and sign(//) = ±1 with a 
flat distribution. Using Micromegas we compute all observables for each sampling point and 
select those rare cases that reproduce all experimental data within 5a errors. 

We do not use any special technique such as Markov chains, which makes the scanning more 
efficient, but that might lead to missing some local best-fit regions. Using the computing power 
of the Baltic Grid we get about 200, 000 good points that satisfy all the experimental criteria. 
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Figure 4: The (M-qui^si) plane in the CMSSM. In the left panel we show the global fit: the 
yellow regions surrounded by continuous contours are the best fit including the Xenon 1 00 and 
LHC data, at 68,95,99.7% confidence levels for 2 d.o.f. The red (blue) regions surrounded by 
dashed contours are the corresponding regions now excluded by Xenon 100 (LHC). In the right 
panel we show points with A% 2 < 4 2 , colored according to the DM annihilation mechanism. 
The red dots in the upper region excluded by the Xenon 100 correspond to the "well-tempered" 
neutralino, green via the heavy Higgs resonance, cyan via neutral Higgses with tan -enhanced 
couplings, blue via slepton co- annihilations, magenta via stop co- annihilations. 

For such points we compute the global \ 2 using the data summarized in Table 2: 

X — Xsm parameters + Xobservables + XlHC + XXenonlOO ■ 

(14) 

Notice that we do not include electroweak precision observables, that do not have a significant 
impact on the result, and that are not well approximated simply by the oblique parameters 
[53]. 

In this work we perform a purely phenomeno logical fit, considering values of the sparti- 
cles masses significantly above the weak scale (up to 4TeV), ignoring the theoretical issue of 
naturalness (see [54, 55] for a recent analysis). Technically, this is achieved as follows: when 
plotting the \ 2 as function of one or two parameters, we minimize it with respect to all other 
parameters. The fit is mainly driven by the DM abundance and by the apparent anomaly in 
the anomalous magnetic moment of the muon, and agrees with the fits in [40, 41, 56]. Given 
that it might not be a real anomaly, we also show regions at relatively high confidence levels. 

We keep the nuclear matrix elements and the DM local density fixed to their default values 
in Micromegas, as already discussed in the previous section. 

5.2 The CMSSM results 

Fig. 4a shows our global CMSSM fit for the DM mass M DM and spin- independent DM-nucleus 
cross section <7si measured by XenonIOO experiment. The yellow regions surrounded by con- 
tinuous contours are the best fit regions including the XenonIOO and LHC data, at 1,2 and 
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Figure 5: Global CMSSM fits in (m , M 1 / 2 ) (left) and in (//, M 1 / 2 ) (right) planes. The red (blue) 
contours are excluded by Xenon 100 (LHC), other details are as in Fig. \a. 



3a level (68,95,99.7% confidence levels for 2 d.o.f). We also show, as red regions surrounded 
by dashed contours, the previous best-fit regions at the same confidence levels now excluded 
by XenonIOO at more than 3a. Obviously, such excluded regions lie around the XenonIOO 
exclusion bound at 90% confidence level (the continuous curve in the figure). 

Within the CMSSM, thermal freeze-out of neutralino DM can reproduce the observed DM 
cosmological abundance according to a few qualitatively distinct mechanisms, that correspond 
to different fine-tunings. To interpret this result we therefore discriminate such distinct cases, 
plotting in Fig. 4b the points of the CMSSM parameter space (also imposing a reasonably good 
global fit, Ax 2 < 4 2 ) colored according to their dominant DM annihilation mechanisms. We 
identify the following distinctive regions: 

1. Red points ("arm" in the upper part of the plot) have ~ Mi such that the lightest 
neutralino has a significant Higgsino component ("well-tempered" neutralino [20]). In the 
CMSSM this happens in the fine-tuned region with large mo and small fi in the so-called 
"focus-point" region [21]. The double structure of the "arm" in Fig. 4b shows that both 
signs of the fi parameter can give the correct DM relic abundance while /i > is favored 
by the (g — 2) M anomaly. Comparison with Fig. 4a shows that the new XenonIOO data 
excludes or disfavors this scenario. 

2. Blue points (around the center of the plot) have M dm , corresponding to slepton 
co-annihilations. They provide the best global fit. 

3. Green points (around the bottom-right part of the figure) have Mdm ~ corre- 
sponding to the resonant DM annihilations mediated by the heavy CP-odd Higgs boson 
A. The light-Higgs resonance (that arises for Mdm ~ rnh/2 and consequently implies 
a relatively light gluino) has now been excluded by the new bounds on the gluino mass 
obtained by the CMS and ATLAS collaborations after 1.1 /fb of data. 
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Figure 6: Global CMSSM fits in (Mg,o~si) (left) and in (m^asi) (right) planes. The red (blue) 
contours are excluded by Xenon 100 (LHC), other details are as in Fig. l±a. 

4. Cyan points represent parameters in which one moves away from Higgs resonances (we 
arbitrarily assume a 10% off-degeneracy to discriminate from the previous case) but the 
DM annihilation processes are still mediated via all neutral Higgses in s-channel. A large 
enough Higgs couplings to the SM fermions and to DM is obtained thanks to tan /3 > 45 
and to a non-negligible Higgsino component of the DM neutralino. Consequently some 
of these points have a large direct detection cross section. 

5. Magenta points (at very low ersi) correspond to stop co-annihilations. The latter region 
is/will be further constrained by the LHC data that presently excludes light stops. 

We remark that, within the CMSSM, the LEP and LHC bounds have excluded one more 
mechanism that was considered more plausible because its non-fine-tuned nature: the bino 
annihilations via light slepton exchange in t-channel. In the light of our results the experimental 
data favours slepton co- annihilation as the mechanism for generating neutralino thermal relic 
abundance. This result is a consequence of the fact that the fit is largely dominated by the 
measurement of muon anomalous magnetic moment that requires relatively light sleptons. 

In the light of these considerations, we can now discuss the global CMSSM fits for other 
parameters. Fig. 5 shows some 2-dimensional fits as indicated in the figure. We see that 
XenonIOO strongly disfavors regions with small \\x\ and partly also with small M\/2 (the very 
small allowed region at 3a level for small M\ii below the red excluded region corresponds to the 
Higgs resonance DM annihilation mechanism). Within the CMSSM a small can be realized 
in the "focus-point" scenario that requires a multi-TeV m®. Such a region was previously not 
favored by global fits and is now further disfavored by the XenonIOO bound that disfavors a 
neutralino DM with a significant higgsino component. The fact that small values of below 
400 GeV are now excluded by XenonIOO data is perhaps one of the most important results of 
our work. 

Fig. 6a compares the LHC prospects (here characterized by the approximate gluino mass 
reach), with the prospects of direct detection experiments (characterized by ctsiMdm/Mz, which 
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Figure 7: Global CMS SM fit before (dashed curve) and after (continuous curve) the Xenon 100 
and LHC results. In the upper panels we fit the DM mass (left) the mo (middle left), M\i% 
(middle right), fj, (right) parameters; in the lower panels we fit tan/3 (left), the Higgs mass 
(middle left), the muon g — 2 (middle right) and the fine-tuning parameter (right). 

is the only parameter probed by experiments when Mqm is much heavier than the target 
nucleus). We see that LHC excluded the Higgs resonance region, that within the CMSSM 
necessarily has a light gluino with mass M3 f» 400 GeV, and reached the best-fit region. On 
the other hand, Fig. 6b shows that, after XenonIOO, there is a clear correlation between the 
spin- independent DM direct detection cross section and Although the XenonIOO results 
exclude only a very small region of the parameter space at low values of m^, this implies that 
the CMSSM charged Higgs boson mass must exceed M H + > 400 GeV and is inaccessible at the 
7 TeV LHC. 

Finally, Fig. 7 shows our fits for the DM mass and for a few key parameters: mo, Mi/ 2 , tan f3, 
rrih and 2/i 2 /M§ (which is a simple measure of fine-tuning, closely related to the fi term). In all 
cases new data had a significant impact: the exclusion of lighter sparticles, moves the best fit 
to higher sparticle masses and consequently to higher tan (3 in order to fit the g — 2 anomaly. As 
shown in the g — 2 panel, there starts to be some tension between its experimental measurement 
(error bands) and the CMSSM predictions. Adding the LHC and XenonIOO bounds, the global 
X 2 at the best-fit point worsens by A% 2 m 4.3. Furthermore, the lower right panel shows that 
the fine-tuning grows up to the few hundred range (for a dedicated discussion see the updated 
version of [54]). 

5.3 "Well-tempered" neutralino 

As we have shown, within the CMSSM, the XenonIOO result has interesting implications on 
the "well-tempered" neutralino [20]. Here we study this scenario of DM generation in SUSY 
models in a model-independent way. Within the MSSM, the neutralino is a mixed state of 
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well tempered bino/higgsino, tan ft = 3 



well tempered bino/higgsino, tan = 10 




Figure 8: Model independent "well-tempered" neutralino scenario. The 3a range for the cos- 
mological DM abundance is reproduced within the green strip. The gray region is excluded by 
Xenon 100 [16]. 

bino, wino, higgsino. None of them, in a pure state, allows thermal DM with a weak scale mass 
M DM « M z : 

• The pure higgsino couples to the Z too strongly, such that for Mdm ~ Mz its ther- 
mal abundance is too low (the cosmological abundance is obtained for a heavy M DM m 
1 TeV higgsino); furthermore, for the same reason it is experimentally excluded by direct 
searches. 

• The pure wino similarly has too much co-annihilations with charged winos, such that for 
Mdm ~ Mz its thermal abundance is too low (the cosmological abundance is obtained 
for a heavier Mdm ~ 2.7 TeV taking into account electroweak Sommerfeld effects [57]). 
Contrary to the previous case, having no coupling to the Z it is allowed by direct searches. 

• The pure bino, instead, has no couplings and no co-annihilations, such that its cosmolog- 
ical abundance would be too high. 

Given that the bino has opposite problems with respect to the higgsino or the wino, it is 
possible to find a good DM candidate by appropriately mixing them [20]. A mixed bino/wino 
still has no couplings to the Z, such that it is not interesting for direct detection; furthermore it 
requires M\ M2 at the weak scale and is not compatible with unification of gaugino masses, 
Mi « M 2 « M 3 at the GUT scale. 

We thereby focus on a mixed bino/higgsino. In the limit where we can ignore all other 
heavier sparticles, its phenomenology is fully described by 3 parameters: the bino mass term 
Mi, the higgsino mass term fi (we assume them to be positive) and tan (3. The observed thermal 
relic DM abundance is reproduced in the green strip in Fig. 8 (left panel for tan = 3 and right 
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Figure 9: The pMSSM fit (left panel) in (Mdmj csi) plane and the corresponding DM generation 
mechanisms (right panel) in the colour code analogous to Fig. 4 except for three new regions: 
the pink dots denote "well-tempered" bino/wino; the dark red dots the Z -resonance at Mdm ~ 
M z /2; the orange dots denote the light-higgs resonance at M DM « 

panel for tan f3 = 10). The region with Mi « fi m Mz was allowed, but its large direct detection 
cross section is now disfavored by XenonIOO (gray region). An improvement of the XenonIOO 
bound by a factor of few would fully exclude the whole "well-tempered" neutralino scenario, 
unless the local DM density or the nuclear matrix element / of eq. (12) are significantly lower 
than what is assumed in our computation. 

The minor tilt at Mdm ~ mt is due to the top quark threshold. At lower masses, the cosmo- 
logically allowed region of Fig. 8 is affected by the Z and Higgs resonances (2M DM = M z or m h 
respectively, indicated as red curves). At larger masses, the "well-tempered" neutralino region 
terminates at fi ~ 1 TeV, where the (almost) pure higgsino becomes a good DM candidate. 

5.4 Supersymmetry beyond the CMSSM: pMSSM 

In order to relax the theoretical constraints of the CMSSM on the SUSY mass spectrum, we also 
performed a similar analysis for free low energy SUSY parameters of the MSSM. The low energy 
phenomenological MSSM, pMSSM, is characterized by the three gaugino masses Mi,M2, Ms, 
the higgsino mass /x, the common squark masses rriq, the left and right-handed slepton masses 
m^,m^, the Higgs mass parameter and by tan/3. We randomly scan all these parameters. 

Fig. 9 shows our results, in the same notations as the corresponding CMSSM figure, Fig. 4. 
Because the squark and slepton masses are not related any more via the GUT relations, light 
sleptons become available for a good fit of (g — 2) M that dominates the fit. At the same 
time, heavy squarks can suppress any new contributions to b — > sj. Therefore a wider range 
of parameters becomes allowed by the global fit, with still a preference for relatively light 
DM mass, that again mostly corresponds to neutralino/slepton coannihilations (blue dots in 
Fig. 9b). Again the region disfavored by the new XenonIOO data mostly corresponds to the 
"well-tempered" bino/higgsino (upper red dots in Fig. 9b). Two new regions appear: "well- 
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tempered" bino/wino (pink dots, at any mass with <tsi below the XenonIOO bound), and 
Z- resonance (dark red dots at There are now significant overlaps among the 

various kinds of DM neutralinos and the mechanisms for generating the DM relic abundance 
are not well separated. Although the pMSSM parameter space is more complicated than the 
CMSSM one, the results are qualitatively similar. 

6 Conclusions 

We have performed a fit to the new XenonIOO data and, based on that, analyzed several 
well motivated DM scenarios. We show that XenonIOO disfavors the Inelastic Dark Matter 
interpretation and the light Dark Matter interpretation of the DAMA/LIBRA claim as well as 
other hints for light DM by CoGeNT and CDMS II. 

The first results of the XenonIOO experiment exclude part of the parameter space of Dark 
Matter models coupled to the Higgs boson. The constrained version of the scalar singlet dark 
matter model (where the DM mass is predicted in terms of DM coupling to the Higgs boson), 
previously favored by the CDMS hint, gets disfavored. Such model now needs an independent 
mass term for the scalar singlet around the weak scale, with its associated naturalness problem. 

In the context of supersymmetry, we show model independently that the XenonIOO data 
disfavors the generic scenario of "well-tempered" bino/higgsino proposed as a way to get neu- 
tralino DM around the weak scale. Within the CMSSM, the XenonIOO data excludes small 
small tua regions of the parameter space previously allowed in the global fit at 3cr level. 
As a result the "focus-point" region is ruled out and, most likely, the CMSSM charged Higgs 
boson mass is not accessible at 7 TeV LHC. Furthermore the "Higgs resonance" mechanism for 
neutralino annihilation, which implies Mdm ~ mh/2 and consequently a relatively light gluino, 
has now been probed and excluded by LHC data taken during 2011. 

The remaining CMSSM best fit parameter space at 3a level is relatively compact in which 
the DM relic abundance is generated by the slepton co-annihilation processes. In view of the 
new LHC bounds, the CMSSM best fit region moves to larger DM mass and consequently to 
lower direct detection cross section. The remaining CMSSM parameter space is more fine tuned 
than before. 

We obtained qualitatively similar constraints also for the pMSSM, where there is no con- 
nection between masses of colored and uncolored sparticles. As a consequence the light Higgs 
resonance is still allowed. 

Note added (27/7/2011). The paper has been fully updated at the light of the new results 
presented at the EPS-HEP-2011 conference (21-27 july 2011), including data from ATLAS and CMS 
on supersymmetric particles and about MSSM Higgs bosons, new data from DO, CMS and LHCb about 
B s —7- and new data from Tevatron on the top quark mass. All those results are available in the 

web site http: //eps-hep2011 . eu. The new figure Fig. 10 shows the impact of present XenonIOO and 
LHC data on CMSSM global fits. The allowed CMSSM parameter space has moved to considerably 
large values of the mass parameters thq and M\/2- As a result, the Higgs resonance region of DM 
freeze-out is now completely ruled out. The allowed values of tan /3 are constrained both from below 
(due to larger sparticle masses, larger tan (5 is needed to fit the muon g — 2 anomaly) and from above 
(due to the new constraint on BR(S S — >• fifi), mainly from LHCb). The DM spin independent direct 
detection cross section is now predicted to be a factor of few below the present XenonIOO bound. 
As seen in Fig. 7, the overall fit to the CMSSM parameters has become worse and fine tuning of the 
CMSSM parameters has further increased. 
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Figure 10: Impact of the latest Xenon 100, Tevatron and LHC data with 1.1/fb on CMSSM 
global fits. In yellow the present global fit at 68 and 95% C.L. for 2 dof; in gray the previous 
result without such data. Fig. 10a shows that LHC bounds move the best fit region to higher 
sparticle masses; then a higher tan/3 is needed to fit the (g — 2) M (see Fig. 10b), and a lower 
direct detection cross section is obtained (see Fig. 10c). The slepton co- annihilation mechanism 
of DM freeze-out is favored; while now both the "well-tempered" bino/higgsino mechanism as 
well as the "Higgs resonance" mechanism are no longer allowed by the global fit. 

Note added (18/4/2011). After our fits to the new XenonIOO data were performed and studies 
of iDM completed, a similar dedicated study appeared by the Xenon Collaboration itself [58]. Our 
results agree in that the iDM is strongly disfavoured as an explanation to the DAMA/LIBRA annual 
modulation signal. 
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Figure 11: Update of fig. 3. Predictions of the scalar singlet model. 

Addendum: 2012 data 

We update our previous results including the new data released during July 2012 by the 
XenonIOO collaboration [1] (sensitivity to asi improved by a factor of 5, see fig. 11) and by the 
ATLAS and CMS collaborations (stronger bounds on sparticle masses, and the measurement 
of the Higgs mass: m h = 125.5 ±0.5 GeV [2]). 

Fig. 11 shows that the scalar singlet DM model considered in section 4 now survives only if 
DM is heavier than about 100 GeV. In particular this excludes the sub-case where all the DM 
mass comes from the Higgs vacuum expectation value (green dot in fig. 11). 

Coming to supersymmetry, fig. 12 shows that the well-tempered bino/higgsino considered 
in section 5.3 has been excluded. We recomputed the DM detect detection cross section ct S t 
now assuming the measured Higgs mass, and we find that the 'well-tempered' region along 
Mi ~ (green strip in the figure, where the DM thermal relic density matches the observed 
cosmo logical DM density) is excluded. The allowed region at ~ 1 TeV corresponds to the 
pure-Higgsino limit. The allowed region at Mdm ~ ^h/2 corresponds to the Higgs resonance 
region; this possibility is however excluded, in the context of models with gaugino unification, 
because it implies a too light gluino, excluded by LHC bounds. 

Fig. 13 shows in the (Mdm, c"si) plane the CMSSM points that lead to the observed DM 
abundance and are compatible with present data, including the recent Higgs mass measure- 
ment. From a theoretical point of view the CMSSM is a sort of 'spherical cow'; most of its 
parameter space has now been excluded, and we are here focusing on its remaining peculiar 
tails and pieces. In particular, the DM abundance can be explained by alternative fine-tuned 
mechanisms, here plotted as different colours. In decreasing order of <7si, red dots correspond 
to the 'well-tempered' bino/higgsino (excluded), extending into the pure-Higgsino limit. Cyan 
dots correspond to DM annihilations via couplings enhanced by tan j3 ~ 50. Green dots cor- 
respond to annihilations via heavy Higgs resonances with large tan/3. Blue points correspond 
to slepton co-annihilation and magenta points to stop co-annihilation. No real best-fit region 
emerges. 

References [2] ATLAS Collaboration, arXiv:1207.7214. CMS 

Collaboration, arXiv:1207.7235. 

[1] XENONIOO Collaboration, arXiv:1207.5988. 
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well tempered bino/higgsino, tan /3 = 3 



well tempered bino/higgsino, tan /? = 10 




Figure 12: (Update of fig. 8). Model independent "well-tempered" neutralino scenario for = 
125 GeV. The 3cr range for the cosmological DM abundance is reproduced within the green strip. 
The gray region is excluded by Xenon 100 [16]. 
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Figure 13: (Update of fig. 5). The (Mt>m,o~si) plane in the CMSSM. Points with Ax 2 < 5 2 ; 
colored according to the DM annihilation mechanism. 
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